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cn . Abstract 

O. 

^ . in and around a cluster. We show that the evolution of the disk surface density 

affects the efficiency of ram-pressure stripping of galaxies at z^l. We also consider 
Qh' the saturation of thermal conduction in detail and show that it cannot be ignored at 

O ' larger radii of a cluster, which makes the time-scale of the evaporation larger. Both 

'^ \ the ram-pressure stripping and evaporation could affect the evolution of galaxies even 

around a cluster. In particular, the observed gradual decline of the star-formation 



'i 



In this paper, we focus on ram-pressure stripping and evaporation of disk galaxies 



rates of galaxies in and around clusters could be explained by evaporation without 
resorting to speculative strangulation (stripping of warm gas in galactic halos). 

Key words: galaxies: clusters: general — galaxies: evolution — galaxies: high- 
redshift — galaxies: interactions 



1. Introduction 

Clusters of galaxies in the redshift range of 0.2-0.5 often exhibit an overabundance, 
relative to present-day clusters, of blue galaxies (Butcher, Oemler 1978, 1984). This star- 
formation activity is often called the Butcher-Oemler effect. Subsequent studies have confirmed 
this trend (Couch, Sharpies 1987; Rakos, Schombert 1995; Lubin 1996; Margoniner, de Carvalho 
2000; Ellingson et al. 2001; Goto et al. 2003a). 
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On the other hand, a correlation has been known between galaxy morphology and the 
local environment. Dressier (1980) studied 55 nearby galaxy clusters and found that fractions 
of early-type galaxies increase and those of late-type galaxies decrease with increasing local 
galaxy density in the clusters. Subsequent studies have confirmed the morphology-density 
relation, that is, early-type galaxies are dominated in inner region of clusters where the density 
is high and the fraction decreases toward the outside of the clusters (Whitmore, Gilmore 1991; 
Whitmore et al. 1993). 

Recently, we have gradually understood that the above phenomena are related to each 
other. Dressier et al. (1994) and Couch et al. (1998) found that most of the blue galaxies 
observed as the Butcher-Oemler effect are normal spirals with active star formation. Dressier 
et al. (1997) studied 10 clusters at z ~ 0.5, and found the morphology-density relation at these 
redshifts. However, they also found that SO fractions are much smaller than those in nearby 
clusters. The low fractions of SO galaxies have also been observed by others (Fasano et al. 2000). 
In many clusters, their galaxy population gradually changes from a red, evolved, early-type 
population in the inner part of the clusters to a progressively blue, later-type population in 
the extensive outer envelope of the clusters (Abraham et al. 1996; Balogh et al. 1997; Rakos et 
al. 1997; Oemler et al. 1997; Small et al. 1998; Couch et al. 1998; van Dokkum et al. 1998; Goto 
et al. 2004). These observations suggest that the blue, normal spirals observed in high-redshift 
clusters were originally field galaxies; they fell into clusters and evolved into the non-blue SO 
galaxies observed in nearby clusters. 

Several mechanisms have been proposed that can lead to color and morphologi- 
cal transformations between galaxy classes in clusters, such as, galaxy mergers (Toomre, 
Toomre 1972), tides from the cluster potential (Byrd, Valtonen 1990; Fujita 1998), tidal 
interactions between galaxies (Moore et al. 1996), ram-pressure stripping (Gunn, Gott 
1972; Takeda et al. 1984; Gaetz et al. 1987; Portnoy et al. 1993; Balsara et al. 1994; Fujita, 
Nagashima 1999; Fujita et al. 1999; Abadi et al. 1999; Mori, Burkert 2000; Fujita 2001a), evap- 
oration (Cowie, Songaila 1977), and a gradual decline in the star-formation rate of a galaxy 
owing to the stripping of halo gas (strangulation or suffocation; Larson et al. 1980; Kodama et 
al. 2001; Bekki et al. 2002). 

In this paper, we focus on ram-pressure stripping and evaporation following Fujita (2004, 
hereafter Paper I). We mostly consider the star-formation history of galaxies, and do not treat 
the morphological transition in detail. We consider the redshift evolution of the disk-size and the 
surface density, which were not considered in Paper I. Moreover, we consider the evaporation 
in detail, while paying attention to the saturation of thermal conduction. We mainly treat 
the environmental effects on galaxies during the first infall into a cluster; we do not consider 
their long-term evolution. Thus, the high-redshift galaxies we investigate may not be direct 
progenitors of galaxies at 2; ~ 0. Although we often use the word 'evolution' from now on, it 
does not mean the evolution of 'a particular galaxy'. Instead, we discuss the differences of the 



average properties of galaxies at low and high redshifts. Although the models presented in this 
paper could be included in complicated semi-analytic models of galaxy formation, it would be 
instructive to study the characteristics of the environmental effects using simple models before 
we consider such a semi-analytic approach. 

This paper is organized as follows. In section 2, we summarize our models. In section 3 
we give the results of our calculations, and compare them with observations in section 4. 
Conclusions are given in section 5. As a cosmological model, we consider a cold dark-matter 
model with a non-zero cosmological constant (ACDM model). The cosmological parameters are 
h = 0.7, where the Hubble constant is given by Hq = lOOh km s~^ Mpc~^, Qq = 0.25, Aq = 0.75, 
and as = 0.8. 

2. Models 

2.1. The Growth of Clusters 

The typical mass of progenitors of a cluster can be derived from the extended Press- 
Schechter model (EPS) (Bower 1991; Bond et al. 1991; Lacey, Cole 1993) and its further 
extension (Fujita et al. 2002). The latter is a Press-Schechter model including the effect of 
spatial correlations among initial density fluctuations (SPS). The models are summarized in 
Paper I. 

Cluster progenitors can be classified into two categories: one is the main cluster and the 
others are subclusters. The main cluster is the progenitor that was located near to the center 
of the current cluster, and had a much larger mass than the other progenitors. Subclusters are 
progenitors other than the main cluster that were located in the vicinity of the main cluster. 
Since the SPS model includes information about the spatial correlation of the initial density 
fluctuations of the universe, we can separately derive the typical masses of the main cluster and 
the subclusters. Moreover, the SPS model includes the EPS model. As shown in Paper I, we 
can define the typical mass of the main cluster at redshift z, which becomes a cluster of mass 
Mq at a later time, zq{< z), as 

J^';;^^^MPEPs{M,z\Mo,zo)dM 
lMLPEPs{M,z\Mo,zo)dM 

where Mmin is the lower cutoff mass, and Peps{M,z\Mo,zo) is the conditional probability based 
on the EPS model that a particle that resides in an object ('halo') of mass Mq at redshift 
Zo is contained in a smaller halo of mass M ~ Mi + dM at redshift z [z > zq). We choose 
Afmin = 10^ Mq, whlch correspouds to the mass of dwarf galaxies. The definition of Peps is 
shown in equation (8) in Paper I. 

The typical mass of the subclusters is given by 

!^l^MPsPs{M,z\Mo,zo)dM 
J^t^PsPsiM,z\Mo,Zo)dM 
3 



Mpfs\z\Mo,zq) — j^^ — - — ^^^ , ^^ — "Tur^ ' y^) 



Msps[z\Mo,Zo) - — ]g^ — ,,^ , ,^ riTT' ' l^-) 



where Psps{M,z\Mq,zq) is the conditional probabihty based on the SPS modeL The definition 
of Psps is shown in equation (7) in Paper I. We define the radius of the region that later becomes 
a cluster of mass Mq aX z = zq as Rq. Following Paper I, we consider the subclusters that were 
initially located at 0.7-Ro <r < R^'m. the precluster region. We refer to the inner radius as Rin- 

2.2. Ram-Pressure Stripping 

We adopt the ram-pressure stripping model of Paper I. In the following sections, we often 
refer to a relatively large dark halo containing galaxies and gas as a 'cluster'. This 'cluster' 
includes the main cluster, the subclusters, and so on. 

We assume that a cluster is spherically symmetric and the density distribution of the 
dark matter is 

Pr^ir) = pmvir/r^i,)'" , (3) 

where pmv and a are constants, rvir is the virial radius of the cluster, and r is the distance from 
the cluster center. We choose a = 2.4 and determine pmv and rvir by a spherical collapse model 
following Paper I. We note that the average mass density of a cluster increases toward high 
redshift. For example, it is proportional to (1 + zY for the Einstein-de Sitter universe. 

We ignore the self-gravity of the ICM, and consider two ICM mass distributions. When 
the ICM is not heated by anything other than the gravity of the cluster, the distribution is 
written as 

where rc/rvir = 0.1. We call this model 'the non-heated ICM model'. In this model, we assume 
that the ICM temperature equals to the virial temperature (Ticm = ^vir) • 

However, at least for nearby clusters and groups, an entropy excess of the ICM (an 
entropy fioor) has been observed in X-rays in the central regions. This indicates that the ICM 
has been heated by some sources, such as AGNs or supernovae, in addition to the gravity of the 
clusters and groups (Ponman et al. 1999; Lloyd-Davies et al. 2000; Mulchaey 2000; Mulchaey 
et al. 2003). Thus, we construct 'the heated ICM model'. We assume that the ICM is heated 
before cluster formation. Although there is a debate about whether the heating took place 
before or after cluster formation (Fujita 2001b; Yamada, Fujita 2001; Babul et al. 2002; Voit 
et al. 2002), the following results would not be much different, even if the ICM is heated 
after cluster formation (Loewenstein 2000). If the ICM is heated non-gravitationally before 
cluster formation, the final distribution depends on the virial temperature of the cluster, T^i^. 
If ^vir > ^0; "wc assumc that a shock forms near the virial radius of the cluster. In this study, 
we assume Tq = 0.8 keV from X-ray observations (Fujita, Takahara 2000; Paper I). The ICM 
distribution is given by 

, , [l + (r/re)^]'36/2 

PiCM(r) - PiCM,vir|^^^^-y^r)ip7^ ' (5) 
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where b = (2.4/3)Tvir/TicM and Tjcm is the temperature of the ICM, which is given by Tjcm = 
Tvir + T}). If Tvir < To keV, a shock does not form, but the gas accreted by a cluster adiabatically 
falls into the cluster. The ICM density and temperature profiles are approximately given by 



Picuir) — picM,^ 



A \ r 



3/2 



4 
Ticuir) = — 



A \ r 



(6) 

(7) 



where A is the constant determined by the adiabat of the ICM (Balogh et al. 1999a; Paper I). 
Since equations (6) and (7) diverge at r = 0, we take their values at r = O.lrvir as the central 
values. In equations (4), (5), and (6), the normalizations of the ICM profile are given by the 
observed ICM fraction of clusters or the rate of gas accretion to clusters (Paper I). If the non- 
gravitational heating makes the accretion time larger than the lifetime of a cluster, the cluster 
cannot accrete much gas and the gas fraction is smaller than the average in the universe. This 
happens at high redshifts {z ^ 1-2) in our models. 

We consider a radially infalling disk galaxy from the turnaround radius of a cluster 
(2rvir)- As the velocity of the galaxy increases, the ram- pressure from the ICM also increases. 
The condition of ram-pressure stripping is 



PlCM^'rel > 27rG'S^SHI 

= 2.1 X lO^^Myncm^^ 



"^rot 



220 km s' 



VlOkpc; V8x 1020mHcm-2y ' ^^ 

where Urei is the relative velocity between the galaxy and the ICM, E* is the gravitational 
surface mass density, Shi is the surface density of the H I gas, f rot is the rotation velocity, and 
Tgai is the characteristic radius of the galaxy (Gunn, Gott 1972; Fujita, Nagashima 1999). We 
define the cluster radius at which condition (8) is satisfied for the first time as the stripping 
radius, rgf Since we assume that the ICM is nearly in pressure equilibrium for r < Tvir, the 
relative velocity, frei, is equivalent to the velocity of the galaxy relative to the cluster, v, for 
r < r^ij.. 

2. 3. Evaporation 

The time-scale of the evaporation of cold gas in a galaxy is written as 

2 /imn \j^\ 
where ks is the Boltzmann constant, Ticm is the ICM temperature, /i (=0.6) is the mean 
molecular weight, mn is the hydrogen mass, Mcoia is the mass of cold gas in the galaxy, and L 
is the energy flux from the hot ICM surrounding the galaxy via thermal conduction (Paper I). 



We define neutral and molecular gas confined in a galactic disk as cold gas; we do not consider 
the cold gas in a galactic bulge. If the electron mean-free path is smaller than the spatial scale 
of the temperature gradient around the galaxy, the thermal conduction is not saturated and 
the energy flux is given by 

I ^nsat I = 47rr2^iKo7'ic^M > (10) 

where rgai is the galaxy radius, and kq = 5 x 10"'^ erg cm~^ s~^ K~^'^ (Paper I). If the elec- 
tron mean-free path is larger than the spatial scale of the temperature gradient, the thermal 
conduction is saturated and the energy flux is given by 

li^satl = 4vrr|,, X OAn^k^T.cu f ^Ml^") ' ' , (11) 

where ric is the electron number density, and rric is the electron mass (Cowie, McKee 1977). For 
convenience, we often refer to tcond given by L = Lnsat as tcond.nsat and that given by L = Lsat as 
Wd.sat- An actual energy flux is given by |L| = min(|Lnsat|, l^satD- 

2.4- Evolution of Disk Properties 

The most important difference between this study and Paper I is that we consider the 
redshift evolution of disk properties in this paper. We adopted a simple model, discussed in 
Mo et al. (1998), for the galactic disk. 

For a given rotation velocity, the galaxy radius at redshift z is given by 

r,^i{z)=r,,i,o{H[z]/Hor\ (12) 

where rgai,o is the galaxy radius at z = and 

H{z) = Ho[Xo + (1 - Ao - l)o)(l + z)' + 1^0(1 + ^)T^' (13) 

is the Hubble constant at redshift z. The surface density and mass of the disk at redshift z are 
given by 

E,{z) = E,^oH{z)/Ho, (14) 

Mdisk(^) = M^,,^^o{H[z]/Ho)-' , (15) 

where S^q and Mdisk.o are the surface density and mass at z = 0, respectively. As shown in 
figure 1 in Mo et al. (1998), H{z = 1) ~ 2Hq. Thus, the disk radius (surface density) at 2 ~ 1 
is a factor of two smaller (larger) than that at z = 0. 

We give the column density of the Hi gas in the disk and the mass of the cold gas 
before the galaxy is affected by environmental effects as follows. We assume that Shi oc S* for 
simplicity; in the future study, we will not assume this by using a semi-analytic model of galaxy 
formation. Thus, the column density of the H I gas at redshift z is given by 

Eni{z) = Em,oHiz)/Ho, (16) 

where Ehi,o is the column density at z = 0. We also assume that Mcoid oc Mdisk and 
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Fig. 1. (a) Mass evolution of tlic main cluster (dotted line), the typical subcluster (solid line) of the 
LCL, and the group (dashed line), (b) Mass evolution of the main cluster (dotted line) and the typical 
subcluster (solid line) of the HCL. 



Meoid(2) = M,oid,o{H[z]/Ho)-^ , 
where Mcom.o is the mass of the cold gas at z = 0. 

3. Results 
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We considered the evolution of clusters with three masses at z = zq. Two of them are 
the same as those in Paper I. The masses of the two clusters at 2;o = are Mq = 1 x 10^^ Mq 
and 6.7 x 10^'^ Mq. We call the former cluster 'the low-redshift cluster (LCL)', which is studied 
to be compared with clusters observed at 2; ~ 0.5. Actually, at z = 0.5, the mass of the main 
cluster is M^^^ = 3 x 10^"^ Mq, which is close to the masses of well-known clusters observed 
at 2; ~ 0.5 (Schindler 1999). The typical mass of the subclusters of the LCL at z = 0.5 is 
Mvir = 6.7 X 10^3 Mq (figure 1). The latter cluster with Mq = 6.7 x 10^^ M© at ^o = is 
investigated to be compared with the subclusters. Since the mass scale of this cluster is that 
of groups of galaxies, we call this cluster 'the group'. For the group, we considered only the 
evolution of the main cluster. The third cluster that we studied has a mass of Mq = 1 x 10^^ Mq 
at zo = 0.5. This cluster was studied to be compared with recent observations of clusters 
observed at 0.5 ^ 2; ^ 1. We call this cluster 'the high-redshift cluster (HCL)'. 

In figure 1, we present the evolutions of the cluster masses. We expect that when the 
mass of the main cluster satisfies the relation Mvir/Mo = (-Rin/-Ro)^5 the subclusters begin to 
be included in the main cluster and become subhalos. For the parameters we adopted (e.g., 
i?in = 0.7-Ro)5 the subclusters of the LCL and HCL are absorbed by the main clusters, together 
with the galaxies in them, at z ^ 0.4 and 1, respectively. This means that at 2; ~ 0.4 (1) the 
subclusters of the LCL (HCL) are to be observed just outside the main cluster of the LCL 
(HCL). 

Since we are interested in galaxies at high redshift, we investigated relatively large galax- 
ies that can be observed in detail. We consider two model galaxies following Paper I. While 



we fixed the rotation velocities of tlie galaxies t^rot, we changed the radius, surface density, and 
Hi column density according to equations (12), (14), and (16). The parameters for the bigger 
model galaxy, which is similar to the Milky Way, are Wrot = 220kms~^, 'rgai,o = 10 kpc, and 
Shi,o = 8 X lO^'^TTT-H cm~^. The parameters for the smaller model galaxy, which is similar to 
M33, are Wrot = 105kms~^, rgai,o = 5 kpc, and Shi.o = 14 x IO^'^ttt-h cm~^. These are the values 
at z = even for the galaxies in the HCL. From now on, this disk evolution is considered unless 
otherwise mentioned. 

3.1. Ram-Pressure Stripping 

First, we consider the time-scale of ram-pressure stripping. Although equation (8) is the 
condition at the representative radius of a galaxy, numerical simulations performed by Abadi, 
Moore, and Bower (1999) showed that it can also be applied at 'each radius' of the galaxy. 
Assuming that tirot is the constant and Shi fx S^, the Hi column density has the relation 
Shi oc S^ oc v^^^/r oc r~^, where f is the distance from the galaxy center. Thus, equation (8) 
shows that the ram-pressure required for stripping at f = 5 kpc is 4-times larger than that at 
f = 10 kpc. Figure 2 shows the time elapsed since the cold gas at f = 10 kpc is stripped until 
that at f = 5 kpc is stripped, Atio_5, for the radially infalling bigger galaxy {v^-ot = 220kms~^). 
For f ^ 5 kpc, the effect of the galactic bulge would be important; since we are interested in the 
star-formation activities in galactic disks, we ignore the inner regions. For a comparison, the 
rotation time of the galaxy at f = 10 kpc is shown in figure 2. The minimum time-scale of the 
evolution of a galaxy is expected to be the rotation time (trot)- For example, the passage of a 
galactic arm through a molecular cloud could stimulate star formation. Thus, if Atio_5 ^ trot, 
the ram-pressure stripping can be regarded as an instantaneous phenomenon for the galaxy, 
which is the case for our model galaxy (figure 2). Thus, the long-term evolution of the galaxy 
does not need to be considered. Moreover, Atio_5 is much smaller than the crossing time of 
the galaxy in the cluster. If turbulent and viscous stripping is considered, the stripping could 
be even faster (Quilis et al. 2000). Since the picM and Vj-d rapidly decrease as the distance from 
the cluster center increases, the ram-pressure becomes less effective in the outer region of the 
cluster. Thus, at cluster radii r > rgt, where rst is defined by the pressure balance at the galactic 
radius f = rgai [equation (8)], it is unlikely that the ram-pressure directly affects the evolution 
of the galactic disk. Although galactic gas at f > rgai might be affected by ram-pressure at 
r > Tst, the effect should be included in 'strangulation', which is considered in Paper I. 

Since cold gas in a galactic disk is almost instantaneously stripped by ram-pressure, 
^st/''"vir should be related to the fraction of galaxies affected by ram-pressure stripping in a 
cluster. The evolutions of rgt/rvir are shown in figures 3 (the bigger galaxy) and 4 (the smaller 
galaxy). Compared to the non-heated ICM model, rst/rvh decreases faster toward higher red- 
shift in the heated ICM model. This is because the non-gravitational heating reduces the 
ICM density and the ram-pressure on galaxies in the inner part of a cluster. In the heated 
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Fig. 2. Time-scale of ram-pressure stripping, Aiio-5 for the bigger galaxy in the (a) main cluster of the 
LCL, (b) subclustcr of the LCL, (c) group, (d) main cluster of the HCL, and (c) subclustcr of the HCL. 
The solid lines are the results of the non-heated ICM model and the dotted lines are those of the heated 
ICM model. The dashed- line shows the rotation time of the galaxy (irot)- 

ICM models, the changes of the slopes correspond to the transformations of the assumed ICM 
distributions (see Paper I). The differences between figures 3 and 4 are small, which shows 
that the differences of the galaxy properties do not much affect the ram-pressure stripping. Of 
course, for much smaller galaxies, r/rgt should be much different from those in figures 3 and 4. 
However, those small galaxies are difficult to be observed in detail at high redshifts. 

In figures 3 and 4, we also present the evolutions when the disk properties do not evolve, 
which were shown in figures 2 and 3 in Paper I; for z ^ 0.5, when the disk properties evolve, 
?"st/''"vir is not much different from that when they do not evolve. For z^l, however, the former 
is much smaller than the latter. This is because in the former case, the surface density of the 
disk is larger at a higher redshift, and the disk is more robust against ram-pressure stripping 
[equation (14)]. As discussed in Fujita (2001a), the ram-pressure from the ICM is larger at 
higher redshifts, because the average mass density of clusters at higher redshifts is larger than 
that of clusters at lower redshifts. Thus, the results show that the effect of the larger disk 
surface density dominates that of the larger ram-pressure at high redshifts. 

For the HCL, rgt/Tvir is larger than that for the LCL for a given redshift, because the 
mass of the HCL is larger than that of the LCL, and the typical galaxy velocity in the HCL 
is larger than that in the LCL (Paper I). On the other hand, for a given mass, ram-pressure 
stripping is more effective in higher redshift clusters. The mass of the main cluster of the HCL 
at 2; ~ 1 is almost the same as that of the LCL at 2; ~ 0.5 (~ 3 x 10^^ Mq; figure 1). However, 
"Tst/f^viv is larger and ram-pressure stripping is more effective in the former. 

Tormen, Moscardini, and Yoshida (2004) studied the orbital properties of galaxies in 
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Fig. 3. Evolution of the stripping radii, rgt, normalized by the virial radii, ryir, for the bigger galaxy in 
the (a) main cluster of the LCL, (b) subcluster of the LCL, (c) group, (d) main cluster of the HCL, and (e) 
subclustcr of the HCL. Solid lines arc the results of the non- heated ICM model and dotted lines are those 
of the heated ICM model. Bold and thin lines are those when the evolution of a galactic disk is considered 
and those when it is not considered, respectively. For r < rgt, ram-pressure stripping is effective. 
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Fig. 4. Same as figure 3, but for the smaller galaxy. 
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massive clusters at z < 0.8 by numerical simulations. Since they showed that the typical peri- 
centric radius of galaxies for the first passage is about O.Srvir, ram-pressure stripping is effective 
if rst/Tvir ^ 0.3. Here, we assume that the results of Tormen, Moscardini, and Yoshida (2004) 
can be applied to our model clusters and the progenitors. If the disk evolution is considered and 
non-gravitational heating is not considered, ram-pressure stripping is effective at z ^ 1 for the 
main cluster of the LCL, z ^ 2.5 for the main cluster of the HCL, and z ^ 1 for the subclusters 
of the HCL, while it is not effective for the subclusters of LCL and the group regardless of z 
(figures 3 and 4). Figures 3 and 4 also show that if both the disk evolution and non-gravitational 
heating are considered, ram-pressure stripping is effective at 2; ^ 0.5 for the main cluster of the 
LCL, and z ^ 1.5 for the main cluster of the HCL, while it is not effective for the subclusters of 
the LCL and HCL, and the group regardless of z. We note that in Paper I the tidal force from 
the main cluster and the resultant shift of the orbit of a galaxy in the subcluster are estimated; 
the estimated peri-centric radius of the galaxy is ~ 0.2rvir. Thus, the assumed threshold of 
ram-pressure stripping {r^^jr^-„ ^ 0.3) seems to be reasonable even for the subclusters. For the 
subclusters of the LCL, for example, since 0.2 ^ rgt/rvir ^ 0.3 at z ^ 1 in the non-heated ICM 
model , we should rather say that the ram-pressure stripping is marginally effective at 2; ^ 1 in 
the non-heated ICM model. 

3.2. Evaporation 

In Paper I, we mainly investigated the evaporation effect on galaxies in clusters (or the 
progenitors) that have not been heated non-gravitationally. In this study, we also focus on 
the evaporation in clusters that have been heated non-gravitationally. Moreover, we study 
the saturation effect of thermal conduction on the evaporation time-scale in detail. Following 
Paper I, we assume that Mcom.o = 5 x 10^ Mq for the bigger galaxy, and Mcom.o = 4 x 10^ Mq 
for the smaller galaxy. 

In figures 5 and 6, we present tcond,nsat and tcond,sat at r = and rvir for the bigger galaxy; 
figure 5 is for the non-heated ICM model, and figure 6 is for the heated ICM model. Note that 
since tcond.nsat oc Mco\d/fgai for the unsaturated case and tcond,sat oc Mcoid/''"gai for the saturated 
case [equations (9), (10), and (11)], tcond,nsat (^cond,sat) for the smaller galaxy is 1.6 (3.2) times 
larger at a given redshift for our parameters. Thus, the differences between the bigger galaxy 
and the smaller galaxy do not much affect the following results. The actual evaporation time- 
scale of a galaxy at a radius r is given by tcond(^) = max[tcond,nsat('^),^cond,sat(^)]- For example, 
at lower redshifts for the main cluster of the LCL, the saturation cannot be ignored at r = r^i,. 

anQ fccond t-cond, sat- 
in the non-heated ICM model, tcond,nsat increases as z increases because tcond,nsat oc T^^y^ 
and TicM decreases rapidly with the mass of the cluster progenitors [equations (9) and (10)]. In 
the non-heated ICM model, tcond,sat(^ = 0) and tcond,sat(^ = ''^vir) decrease slowly as z increases 
(figure 5). This is because their TicM-dependence is small (oc Tj^j^ ) and the increase of rif, at 
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Fig. 5. Evolution of the conduction times for the non-heated ICM model for the bigger galaxy in the 
(a) main cluster of the LCL, (b) subcluster of the LCL, (c) group, (d) main cluster of the HCL, and (e) 
subcluster of the HCL. The solid lines are tcond.nsat at r = (thin lines) and r — r^i^ (thick lines). The dotted 
lines are icond.sat at r = (thin hues) and r = r^i, (thick lines). Note that tcond,nsat(?' = 0) = icond,nsat(?'vir) 
in this figure. 




Fig. 6. Same as figure 5, but for the heated ICM model. 
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high redshift dominates [equations (9) and (11)]. 

In the heated ICM model, the evolution of the conduction time-scales is more compli- 
cated (figure 6). In general, the conduction time-scales in the heated ICM model are smaller 
than those in the non-heated ICM model because of larger Tjcu owing to non-gravitational 
heating. The detailed behavior of the conduction time-scale when thermal conduction is not 
saturated, tcond.nsat, can be explained as follows. For the main clusters of the LCL and the HCL 
and the subcluster of the HCL, tcond.nsat increases as z increases. This is because Tjcm decreases 
as is the case of the non-heated ICM model. However, the increasing rate is smaller because 
^iCM decreases more slowly owing to non-gravitational heating. For higher redshifts, the viral 
temperature becomes smaller and satisfies the relation Tvir < Tq, where Tq (=0.8 keV) is the 
critical temperature under which the ICM distribution follows an adiabatic accretion model 
(subsection 2.2). For example, T^i^ equals Tq = 0.8 keV a.t z = 1.2 for the main cluster of the 
LCL. Since the cluster is not isothermal in the adiabatic accretion model [equation (7)], tcond,nsat 
bifurcates there. For the subcluster of the LCL and the group, the density and temperature 
profiles follow those predicated by the adiabatic accretion model for z > 0. In these redshift 
ranges, tcond,nsat(r = ^vir) > tcond,nsat(r = 0) bccause Ticuir = Tvir) < Ticuir = 0). In the middle- 
redshift range, 1.2 < z < 1.7 for the main cluster of the LCL for example, the ICM fraction of a 
cluster is the same as that at lower redshift although the cluster is more compact at the higher 
redshifts. Therefore, Ticm increases via adiabatic compression and tcond,nsat ('^ = ''^vir) decreases 
as z increases. At higher redshifts, 2; > 1.7 for the main cluster of the LCL for example, tcond.nsat 
increases as z increases, because the ICM fraction and Ticm decrease. 

The behavior of the saturated conduction time-scale, tcond.sat? at lower redshifts {z ^ 
1.2 for the main cluster of the LCL, for example) can be explained as follows. Since the 
average dark matter and ICM densities increase toward higher redshifts, the increase of tcond,sat 
through the decrease of Ticm is dominated by the decrease of tcond.sat through the increase of n^ 
[equations (9) and (11)] at r = Tvir. However, in our ICM model, for a given gas mass fraction, 
the non-gravitational heating decreases the ICM density in the central region of a cluster. Thus, 
the decrease of tcond,sat is not significant at r = 0. At higher redshifts, once the density and 
temperature profiles follow those predicted by the adiabatic accretion model, the behaviors of 
^cond.sat cau be attributed to the mechanisms that are the same as the case of tcond,nsat. 

In figure 7, we show the evolution of tcond = iiiax(tcond,nsat,^cond,sat) at r = for the bigger 
galaxy in the non-heated ICM model. This figure corresponds to figure 5 in Paper I in which 
the evolution of the disk is not taken into account. As shown in figure 5, tcond,nsat ^ ^cond,sat(^ = 
0) and therefore tcond(^ = 0) = tcond.nsat. The conduction time-scale can be represented by 
^cond oc Mcoid/'^gai [equations (9) and (10)]. Since both M^oid and rgai are proportional to H{z)~^ 
[equations (12) and (17)], tcond is not infiuenced by the effect of the disk evolution. Thus, the 
upper figures are the same as figures 5a-c in Paper I. In figure 7, we also present the Hubble 
time, tn. If ^cond ^ ^h, the cold gas in a galaxy is evaporated soon after the galaxy forms in 
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Fig. 7. Time-scales of tiiermal conduction, tcond(?' = 0) for the non-lieated ICM model for the bigger 
galaxy in the (a) main cluster of the LCL, (b) subcluster of the LCL, (c) group, (d) main cluster of the 
HCL, and (e) subcluster of the HCL. The solid and dashed lines show tcond('' = 0)' = tcow&ij' = 0) +tform 
when Zform = 1 and 2, respectively. The dotted and dot-dashed lines show icond(^ — 0) and the Hubble 
time in, respectively. For t'^^^^^ ^ tn, thermal conduction is effective. 

the cluster or its progenitors. We also show t'^ondi^ = 0) = tcond(^ = 0) +tform, where tform is the 
Hubble time when the galaxy forms ed, z = ^form- If i^cond < ^h at a given redshift, the cold gas 
in the galaxy has been evaporated by the redshift. Note that tcond corresponds to t'^ond when 



Zio 



oo. 



Figure 7 shows that in the main cluster of the LCL (HCL) observed at 2 ~ 0.5 (1), cold 
gas has been evaporated from the galaxies near the cluster center. On the other hand, in the 
subclusters of the LCL (HCL) observed at 2; ~ 0.5 (1), cold gas is evaporating if the galaxies 
formed at z ~ 1-2. The time-scales of the evaporation, tcond, is relatively long (^ 2 Gyr). In 
the group, cold gas is evaporating at 2 ~ 0. 

In figure 8, we present the evolutions of tcond('" = 0) and t'^.^^^{r = 0) for the bigger galaxy 
in the heated ICM model. Since tcond in the heated ICM model is generally smaller than that in 
the non-heated ICM owing to the larger Ticm, ^cond becomes smaller than tu soon after galaxy 
formation. 

A galaxy that first enters a cluster from the outside does not stay at r = r^i,. and the 
orbit is not exactly radial. If we assume that the typical peri-centric radius of galaxies is about 
O.Srvir, ram-pressure stripping is effective for most galaxies if Tst/rvh ^ 0.3 (see subsection 3.1). 
On the other hand, if Tst/r^ij- ^ 0.3, the galaxies would be affected only by evaporation while 
they are orbiting. In this case, the actual time-scale of evaporation, tevap, is given by tcond('" = 
0) < tcvap < ^cond(^ = ''^vir)- Therefore, in the subclusters of the LCL observed at 2 ~ 0.5, the 
bigger galaxy infalling from the outside loses its cold gas with a time-scale of 0.5 < tcvap < 10 Gyr 
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Fig. 8. Same for figure 7, but for the heated ICM model. 

in the heated ICM model (figure 6). In the subclusters of the HCL observed at 2; ~ 1, a galaxy 
loses its cold gas with the time-scale of tgvap ~ 0.5 Gyr in the heated ICM model (figure 6). 

3.3. Comparison between Ram-pressure Stripping and Evaporation 

Although both ram-pressure stripping and evaporation suppress star-formation activities 
in galaxies, they affect the star- format ion activities differently. For position, while the ram- 
pressure stripping is effective only in the central regions of clusters (r ^ Q.br^i^; figures 3 and 4), 
the evaporation is often effective even at r ~ r^i^.. For time-scale, while the ram-pressure 
stripping checks star-formation activities in a very short time (~ 10^ yr; figure 2, see Fujita, 
Nagashima 1999), the evaporation generally affects the star- formation activities more slowly (~ 
10^ yr; figures 5 and 6). The difference in the decline rate of the star-formation activities could 
be discriminated by the spectra of the galaxies. These facts would be useful to observationally 
find whether ram-pressure stripping or evaporation dominates in clusters. In the next section, 
we investigate several specific cases. 

4. Discussion 



In Paper I, the star-formation activities of galaxies in main clusters and the subclusters 
in the vicinity of the main clusters observed at 2; ~ 0.5 are discussed. Since the consideration of 
the disk evolution does not much change rst/rvh for 2; ^ 0.5 for the LCL (figures 3 and 4), the 
conclusions about the effects of ram-pressure stripping on the galaxies do not change either. 
In Paper I, we discussed that the rapid decline of the star-formation rates of galaxies in the 
main cluster of the LCL is not consistent with observations of the CNOC sample of very 
luminous X-ray clusters (Balogh et al. 1999b; Kodama, Bower 2001) and this suggests that 
the star-formation rates have decreased before the galaxies enter the main clusters (Goto et 
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al. 2003b; Goto et al. 2003c). In Paper I, we argued that the 'pre-processing' occurred in the 
subclusters (Zabludoff, Mulchaey 1998; Hashimoto, Oemler 2000; Balogh et ah 2000), which is 
consistent with the observations showing that red galaxies have a clumpy distribution around 
a main-cluster (Kodama et al. 2001). In the subclusters of the LCL at 2; ~ 0.5, ram-pressure 
stripping is marginally effective (0.2 ^ rst/^vir ^ 0.3) in the non-heated ICM, and it is ineffective 
in the heated ICM model (figures 3 and 4). If ram-pressure stripping is the mechanism of the 
pre-processing, most of the cold gas in galaxies may be removed and star-formation activity 
of the galaxies may completely die out before the galaxies enter the main cluster. This may 
be inconsistent with the existence of many blue galaxies in main clusters at 2; ^ 0.5 (Kodama, 
Bower 2001); the observations suggest a slower decline of the star- format ion rates. Thus, the 
heated ICM model is preferable because ram-pressure stripping can be ignored. In the heated 
ICM model, evaporation can be candidates of the mechanism of the pre-processing. The time- 
scale of the evaporation is 0.5 ^ tevap ^ 10 Gyr. Since the time-scale is relatively large, the 
evaporation can be an alternative of strangulation, which is expected to gradually suppress the 
star- formation activities of galaxies (Larson et al. 1980; Kodama et al. 2001; Bekki et al. 2002), 
but is highly speculative (Benson et al. 2000). 

For the HCL observed at 2; ~ 1, ram-pressure stripping is effective in the main cluster 
regardless of non-gravitational heating. As is the case of the LCL at 2 ~ 0.5, the ICM of the 
subclusters must have been heated to avoid ram-pressure stripping in the subclusters. 

At higher redshifts (z ^ 1), the effect of the disk evolution is more significant (figures 3 
and 4), although it is not well-known that disk galaxies with the rotation velocities that we 
studied exist. At these redshifts, most of the subclusters have not been absorbed in the main 
cluster. As long as the ICM has not been heated non-gravitationally, the efficiency of ram- 
pressure stripping decreases only slowly as redshift increases. Thus, the products of ram- 
pressure stripping, such as galaxies that have spectra refiecting rapidly declined star- formation 
rates, may be observed at these redshifts although the number fraction is not as much as 
that at lower redshifts. In the non-heated ICM model, figure 7 shows that the evaporation is 
ineffective (tcond ^ %) for less massive systems such as the subclusters of the LCL (at 2 > 0.7) 
and the HCL (at z > 1.5). Thus, if the evaporation is the mechanism of the pre-processing, 
the star-formation rates of galaxies hardly decline in the vicinity of the clusters at z ^ 1- 
2 except for galaxies in which star- format ion activities decrease rapidly by the ram-pressure 
stripping. On the other hand, if the ICM has been heated, ram-pressure stripping does not 
occur in the high-redshift range (figures 3 and 4). The time-scale of the evaporation is given 
by tcond{r = 0) < tcvap ^ tcond('^vir) and is relatively large; for example, 1 < tevap ^ 4 Gyr at 
z ~ 2 for the subcluster of the LCL (figure 6). Thus, the star-formation rates of galaxies should 
decrease slowly with this time-scale. Thus, the heated and non-heated ICM models may be 
discriminated by future observations of the star-formation history of galaxies in the regions 
surrounding main clusters. However, quantitative predictions based on semi-analytic models 
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or numerical simulations would be required for actual discrimination. We leave this for future 
study. 

5. Conclusions 

We have investigated ram-pressure stripping and evaporation of disk galaxies in and 
around a cluster using analytical models based on a hierarchical clustering scenario. We con- 
sidered the redshift evolution of the size and surface density of the disk. We showed that the 
evolution does not much affect the efficiency of ram-pressure stripping of galaxies for z ^ 0.5, 
but affects for 2 ^ 1. We also considered the saturation of thermal conduction in detail, and 
found that it cannot be ignored at larger radii of a cluster, which makes the time-scale of 
the evaporation larger. Thus, the evaporation has the same effect as strangulation (stripping 
of warm gas in galactic halos) to suppress the star-formation activities of galaxies gradually. 
Observations of galaxies in the vicinity of clusters at z ~ 1 are useful to investigate whether 
non-gravitational heating has occurred or not by 2; ~ 1. 

We are grateful to M. Nagashima and M. Enoki for useful comments. Y. F. was sup- 
ported in part by a Grant-in- Aid from the Ministry of Education, Culture, Sports, Science and 
Technology (14740175). 
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